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SUMMARY 


A pulsed  laser,  remote  crosswind  measurement  system  is 
being  developed  for  remotely  sensing  the  crosswind  in  the 
atmosphere  by  analyzing  the  laser  energy  scattered  from  a 
diffuse  target.  The  system  consists  of  a 100  millijoule, 
double-pulse  Nd:YAG  laser,  a sixteen  element,  pulsed  optical 
receive^and  a processor.  In  order  to  estimate  the  crosswind, 
the  variance,  covariance  and  slope  of  the  time  delayed  co- 
variance  function  for  the  received  intensity  must  be  measured. 
This  has  been  accomplished  by  utilizing  an  optical  receiver 
array  and  a combination  of  space  and  time  averaging  to  esti- 
mate the  statistical  quantities.  Design  and  construction  of 
the  crosswind  measurement  system  has  been  completed  and  it  will 
be  tested  and  optimized  during  the  next  contract  period. 

Theoretical  formulations  for  the  statistics  of  the  received 
intensity  are  needed  for  system  design,  optimization  and  cali- 
bration. During  the  period  covered  by  this  report,  formulations 
were  developed  which  include  the  effects  of  the  log-amplitude 
covariance  as  well  as  the  wave  structure  function  and  represent 
a significant  extension  of  previous  work.  In  addition,  a major 
analytical  breakthrough  was  made  in  our  effort  to  reduce  the 
theoretical  formulations  to  numbers  for  comparison  with  experi- 
mental results,  design  and  calibration.  The  theoretical  formu- 
lations are  being  verified  in  a cw  field  experiment  and  the 
results  thus  far  have  yielded  excellent  agreement  with  this 
theory. 
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I.  INTRODUCTION  AND  BACKGROUND 

In  FY74  Frankford  Arsenal  through  the  Atmospheric 
Sciences  Laboratory  (ECOM)  White  Sands  Missile  Range,  New  Mexico, 
initiated  a program  to  develop  a remote  crosswind  measurement 
system  for  ballistic  applications.  A remote  crosswind  system  is 
needed  that  will  measure  the  average  crosswind  between  gun  and 
target  and  automatically  feed  this  data  to  an  onboard  ballistics 
computer  for  fire  control  correction  of  the  crosswind  error. 

This  error  is  quite  significant,  representing  one  of  the  largest 
remaining  sources  of  combat  vehicle  fire  control  error,  especially 
with  medium  and  slow  muzzle  velocity  ammunitions.* 

The  feasibility  of  remotely  determining  the  trans- 
verse wind  velocity  with  a point-to-point,  cw  laser  system  is 
2 3 A 5-11 

well  established.  ’ ’ In  more  jecent  work,  investigators 

have  considered  the  extension  of  the  technique  to  the  use  of 

multiple  cw  lasers  and  folded  paths  as  well  as  ambient  light 

and  the  use  of  techniques  to  overcome  saturation.  The  technique 

basically  consists  of  relating  the  scintillation  pattern  on  the 

receiver,  to  the  frozen-in,  transverse  motion  of  the  turbulence 

structure  along  the  propagation  path.  The  most  successful  cw 

method  that  has  been  used  involves  the  measurement  of  the  slope 

2 

of  the  log  amplitude  covariance  function  at  zero  time  delay. 

Using  active  cw  sources,  an  acceptable  signal  to  noise 

ratio  can  be  achieved  at  reasonable  transmitter  power  levels. 

However,  the  proposed  application  of  the  work  sponsored  under 

this  project  requires  a single-ended  system  operating  against 

a non-cooperative  target  which  can  be  viewed  as  a perfectly 

12 

diffuse  surface  with  reflectivity  of  approximately  0.1. 

Under  this  constraint,  an  active  cw  system  would  require  ex- 
cessive transmitter  power  levels. 
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The  primary  noise  sources  in  the  system  are  the  Poisson 
fluctuations  of  detector  current  and  thermal  noise  in  the  elec- 
tronics. For  a given  pulse  energy,  as  the  pulse  width  decreases 
the  instantaneous  power  increases.  Consequently,  the  detector 
signal  current  increases  and  eventually  its  Poisson  fluctuation 
becomes  the  dominant  noise  source.  If  the  receiver  is  gated 
synchronously  with  the  transmitted  pulses,  then  it  would 
be  signal  shot  noise  limited.  Therefore,  in  a direct  detection 

optical  system,  the  signal  to  noise  ratio  for  a given  total 

13 

energy  increases  as  the  pulse  duration  is  made  shorter.  Conse- 
quently, a pulsed  system  can  be  made  considerably  more  sensitive 
than  a cw  system. 

In  order  to  employ  a pulsed  laser  and  still  make  a 
temporal  measurement  as  required  for  wind  velocity,  a double- 
pulse technique  has  been  used^’^’1**  whereby  two  Q-switched 
laser  pulses  are  generated  at  a time  interval  on  the  order  of  a 
millisecond.  This  interval  appears  to  be  appropriate  for  the 

2 

employment  of  the  covariance  "slope  at  zero  time  lag"  technique 
and  is  compatible  with  current  double-pulse  laser  technology. 
Initially,  a ruby  laser  was  used  as  a pulsed  source. 

Since  the  pulse  technique  utilizing  a ruby  laser  pre- 
cluded time  averaging  of  the  statistical  quantities  measured, 
spatiai  averaging  had  to  be  employed.  Therefore  t he  equivalent 
of  an  array  of  detectors  was  used  In  order  to  generate  a suf- 
ficient number  of  samples  of  the  instantaneous  scintillation 
pattern  to  use  spatial  averaging.  To  the  extent  that  the 
statistics  of  this  pattern  are  invariant  with  lateral  trans- 
lation at  the  receiver  (a  safe  assumption  with  the  use  of  a 
diffuse  target  as  a secondary  source),  such  spatial  averaging 
will  give  the  desired  result,  and  the  pulsed  technique  is  funda- 
mentally similar  to  the  cw,  time-averaging  approach. 
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Initially  the  work  on  the  pulsed  laser,  remote  cross- 
wind  measurement  system  was  based  on  a technique  and  theory 
developed  for  a cw  system  utilizing  a laser  directly  as  a source. 
This  technique  and  theory,  which  is  based  on  log-amplitude 
statistics,  has  been  verified  experimentally  and  is  valid  for  the 
cw  configuration.  The  basic  cw  technique  is  applicable  to  the 
pulsed  system.  However,  the  detailed  theory  does  not  apply  to 
the  pulsed  configuration.  This  is  not  the  result  of  a differ- 
ence in  the  propagation  characteristics  between  the  pulsed  and 
cw  sources.  It  is  caused  by  the  physical  differences  in  the 
two  systems. 

The  source  for  the  pulsed  system  is  a laser  illuminated, 
incoherent  target  yielding  a source  with  an  undefined  eize  which 
at  least  close  to  the  target  gives  rise  to  a temporally  coherent, 
spatially  incoherent,  Rayleigh  distributed  field  amplitude.  This 
is  in  contrast  to  the  point  source  used  in  the  earlier  active  cw 
systems  and  the  temporally  incoherent  cw  sources  used  in  later 
active  and  passive  systems  vhlch  give  rise  to  log-normal  amplitude 
statist ics. 

Consequently,  the  experimental  work  with  the  ruby 
system  was  concentrated  largely  on  tests  where  the  size  of  the 
source  (at  the  target)  could  be  controlled.16  This  allowed  a 
better  comparison  with  available  theory  by  decoupling  from  the 
experiment  the  effects  of  the  atmosphere  on  the  laser  beam  as  it 
propagates  from  the  transmitter  to  the  target.  The  controlled 
source  size  was  accomplished  by  utilizing  small  disks  of  scotch- 
lite  on  a black  background  as  a target. 

The  results  of  these  tests  indicated  the  need  both  for 
a better  pulsed  source  and  a theory  applicable  to  the  pulsed 
system. 
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A ruby  laser  was  originally  chosen  as  a source  because 
of  cost  considerations  and  the  realization  that  energy  levels  of 
the  order  of  a joule  or  greater  might  be  required.  Subsequent 
analysis  and  experimental  results  indicated  that  energy  levels 
of  the  order  of  100  millijoules  would  yield  an  adequate  signal  to 
noise  ratio  in  the  system.  Consequently,  the  system  was  re- 
designed to  utilize  a Nd:YAG  laser  operatirg  at  1.06  microns. 

It  has  the  following  advantages  over  the  ruby  system:  much 

smaller;  more  rugged;  less  sensitive  to  temperature  changes; 
lighter;  less  power  required;  much  less  RFI;  convection  cooled; 
space-time  averaging  possible. 

Most  of  the  advantages  of  the  Nd:YAG  system  accrue  from 
the  much  lower  threshold  of  the  lasing  medium.  Consequently  it 
requires  a much  lower  Input  energy  per  pulse  (on  the  order  of  10 
In  lieu  of  1,000  joules).  Thin  allows  the  YAG  laser  to  be  double 
pulsed  at  a much  higher  rate  than  the  ruby  laser  and  makes  space- 
time  averaging  possible.  In  addition,  the  lower  Input  energy 
greatly  reduces  radio  frequency  interference  (RKI)  between  the 
laser  and  the  receiver. 

The  use  of  space-time  averaging  allowed  the  size  of 
the  optical  receiver  array  (spatial  averaging)  to  be  reduced 
from  64  elements  (8  x 8)  to  16  elements  (2  x 8).  The  laser 
can  produce  a double  pulse  once  every  100  milliseconds  and  the 
number  of  double  pulses  used  for  a crosswind  measurement  can  be 
adjusted  from  one  to  forty.  Ten  double  pulses  will  normally  be 
used  for  the  time  averaging  which  requires  about  one  second. 

This  yields  160  samples  per  measurement  for  the  YAG  system  as 
compared  to  64  samples  with  the  old  ruby  system.  An  additional 
advantage  of  the  time  averaging  is  that  it  will  tend  to  eliminate 
any  intermittency  of  turbulence  problems.  The  YAG  system  has 
been  constructed  and  will  be  evaluated  during  the  next  contract 
period. 
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Since  a theory  applicable  to  the  pulsed  system  was  not 
available,  a new  theoretical  formulation  for  the  statistics  of 
the  received  intensity  based  on  the  extended  Huygens-Fresnel 
principle  was  developed.  ’ ' This  first  theoretical  effort 
invoked  the  assumption  of  Gaussian  field  statistics  at  the 
receiver.  Because  of  this  assumption,  only  phase  perturbation 
effects,  which  show  up  as  a dependence  on  the  transverse  cor- 
relation scale  (pQ)  appear  in  the  result.  This  was  initially 
thought  to  be  a good  assumption,  at  least  for  low  turbulence 
levels,  on  the  basis  that  the  relation  between  the  two  turbu- 
lence scale  sizes  (pq  and  /L/k)  was  expected  to  be  multiplicative 
and  therefore  the  smaller  scale  - which  is  usually  p in  the 
present  application  - would  be  dominant.  This  approach  has  the 
advantage  of  yielding  a simple  theory,  but  the  disadvantage  of 
requiring  a receiver  to  operate  over  several  decades  of  scale 
size  with  the  scale  size  being  an  exponential  function  of  both 
target  range  and  turbulence  level. 

Our  cw  experiments  have  confirmed  that  at  low  turbu- 
lence levels,  the  Gaussian  assumption  is  indeed  correct  for  the 
variance  and  covariance.  However  at  intermediate  turbulence  levels 
the  scale  size  approaches  /l./k  . This  suggests  that  the  log- 
amplitude  term,  which  is  responsible  for  this  scale  size,  plays 
an  important  part  in  the  statistics  of  the  received  intensity 
and  particularly  with  regard  to  remote  wind  sensing. 

Consequently,  a complete  theoretical  formulation  in- 
cluding the  log-amplitude  term  has  been  developed,  it  shows 
that  the  relationship  between  the  two  scale  sizes  is  additive 
rather  than  multiplicative  as  indicated  by  the  cw  experimental 
data.  This  means  that  /l./k  scale  size  is  dominant  for  our 
application  and  that  the  same  mechanism  which  has  been  used  as 
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a basis  for  the  successful  cw  wind  sensing  work  performed  by 
NOAA  is  available  for  utilization  by  the  pulsed,  diffuse-target 
case.  This  had  not  been  clear  heretofore. 

In  addition,  we  achieved  a major  analytical  breakthrough 
which  greatly  simplified  the  formulations  while  retaining  a 
complete  description  of  the  mechanisms  involved.  We  are  now  able 
to  reduce  the  formulations  to  numbers  for  comparison  with  the  cw 
experimental  results  and  for  design  purposes. 

These  re-'ults  have  very  positive  implications  with  re- 
gard to  the  design  of  a pulsed  remote  wind  sensing  system. 
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II.  1.06  SYSTEM 

During  the  period  of  this  report  a system  was  designed 
and  constructed  to  remotely  sense  the  crosswind  by  utilizing  a 
double  pulsed  neodymium  YAG  laser  as  a transmitter.  The  system 
was  designed  to  take  advantage  of  the  relatively  rapid  repetition 
rate  possible  with  the  Nd:YAG  laser  by  incorporating  temporal 
averaging  in  addition  to  the  spatial  averaging  used  with  the 
ruby  system^* ^ to  obtain  more  accurate  estimates  of  the  sta- 
tistics of  the  received  intensity. 

The  overall  system  is  illustrated  in  Figure  I I— l . It 
consists  of  a double  pulsed  Nd:YAG  laser,  a processor  for  control 
and  data  reduction,  a sixteen  element  optical  receiver,  a four 
channel  analog  to  digital  converter  with  display,  a video  display, 
a wind  and  range  display,  a high  speed  paper  tape  punch,  a tele- 
type and  a control  panel. 

In  order  to  measure  the  crosswind,  the  time  delayed 

statistics  of  the  received  intensity  must  be  measured.  Conse- 

14 

quently,  a double  pulse  laser  must  be  used  as  a source.  The 
laser  used  in  the  system  was  manufactured  by  International  Laser 
Systems  and  Is  basically  their  model  NT-170  modified  for  double 
pulse  operation.  The  specifications  for  the  laser  are  shown  in 
Table  1 I— l and  the  laser  is  shown  in  Figure  I 1-2. 

The  overall  system  control  and  data  processing  neces- 
sary to  carry  out  the  time  averaging  and  perform  on  site  wind 
calculation  was  accomplished  by  using  a Computer  Automation 
LSI-2  minicomputer.  Interface  circuitry  and  software  was  de- 
veloped to  allow  the  minicomputer  to  control  not  only  the  laser 
and  receiver,  but  also  a four  channel  analog  to  digital  con- 
verter and  the  various  input  and  output  dev  chnwn  in  Figure 
II-l. 
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Table  II-l.  Specifications  for  Model  NT-172 
Double-Pulse  Laser  System 


Transmitter: 

Output  Wavelength 
Output  Beam  Divergence 

Beam  Diameter 
Energy  per  Pulse 
Pulsewidth 
Repetition  Rates 
Double-Pulse  Separation 
Operating  Temperature 
Input  Power  Requirement 

Size 

Weight 

Power  Supply : 

Input  Power  Requirement 

Size 

Weight 

Line  Converter: 

Input  Power  Requirement 

Output  Power 

Size 

Weight 


1.064  pm 
0.3  mr 

3.0  in. 

100  mJ 

18  nsec,  nominal 
Single-shot,  1,  5 and  10  pps 
1 millisecond  (adjustable) 
32°  to  120°F 

24  V dc  at  60  A peak,  40  A 
average 

9 x 11  x 25  in. 

43  lb. 


24  V dc 

9 x 11  x 20  in. 
38  lb. 


115  V,  50  to  60  Hz,  10 
24  V dc 

9 x 11  x 20  in. 

70  lb. 


An  operating  system  with  the  capability  to  link  assembly 
level  control  programs  to  Fortran  programs  was  used  to  develop  the 
operating  software.  This  combination  of  machine  language  and 
higher  level  language  allows  the  control  and  data  processing  to  be 
very  flexible.  Scientific  computations  are  carried  out  by  Fortran 
subroutines  which  can  be  written,  modified,  and  debugged  much  more 
readily  than  equivalent  assembly  language  programs.  The  processor, 
wind  and  range  displays,  and  control  panel  are  all  housed  in  a 
common  cabinet  as  shown  in  Figure  I I— 3 . 

The  analog  to  digital  converter  is  designed  to  be  used 
with  propeller  anemometers  and  a Cambell  Scientific  CA-S  space 
averaging  anemometer  to  provide  in  situ  wind  and  turbulence  strength 
data.  This  allows  the  *n  situ  data  to  be  recorded  under  control  of 
the  processor  for  comparison  to  measurements  made  with  the  1.06 
system. 

Electronically,  the  1.06  system  receiver  is  quite  similar 
to  the  ruby  system.^  However,  the  optical  and  detection  portions 
of  the  1.06  system  receiver  are  quite  different  from  that  of  the 
ruby  system. 

The  heart  of  the  receiver  is  a hybrid  array  of  PIN 
silicon  photodiodes.  The  array  consists  of  16  one  millimeter 
diameter  devices  arranged  in  two  horizontal  rows.  The  8 diodes 
in  each  row  are  mounted  on  1.75  mm  centers.  As  shown  in  Figure 
II-4,  a two  lense  optical  system  is  used.  It  is  basically  a 
beam  reducer  that  allows  the  characteristic  scale  sizes  of  the 
turbulence  to  be  reduced  to  match  the  detectors.  The  second 
lense  is  easily  changed  and  provides  a means  of  adjusting  the 
magnification  and  hence  the  effective  array  spacing. 

Results  from  the  cw  experiments  indicate  that  an  array 
spacing  around  0.5/b/k  may  be  optimum.  For  the  initial  experiments 
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Figure  11-4  Receiver  Optical  Design 


at  2000  m (/L/k  = 18.4  mm)  a 60  mm  focal  length  second  lense 
will  be  used.  In  conjunction  with  the  310  mm  objective  lense 
this  produces  a demagnification  of  5.5,  an  effective  detector 
spacing  of  9.625  mm  and  an  effective  detector  aperture  size  of 
5.5  mm.  If  future  experimental  results  indicate  that  adjustment 
of  the  effective  detector  spacing  as  a function  of  range  is 
necessary  to  achieve  the  desirau  accuracy  over  the  full  operating 
range,  then  the  beam  reuucer  could  be  replaced  with  zoom  optics 
that  are  automatically  adjusted  by  the  computer  based  on  the 
measured  range  to  the  target. 

As  shown  in  Figure  11-4,  a field  of  view  limiting 
aperture  is  mounted  at  the  focal  point  of  the  objective  lense. 

This  serves  to  limit  the  amount  of  background  light  incident  on 
the  detectors.  An  aperture  of  1.6  mm  was  chosen  to  limit  the 
lield  of  view  to  5 mill iradians. 

Provision  has  been  made  for  mounting  a removable  beam 
splitter  as  shown  in  Figure  I 1—4 . It  is  a piece  of  uncoated 
crown  glass  at  a 45°  angle  to  the  optical  axis  and  is  designed 
to  divert  about  9%  of  the  light  to  a mean  intensity  detector. 

A mean  intensity  detector  would  provide  a more  accurate  estimate 
of  the  mean  intensity  than  can  be  obtained  from  averaging  the 
16  optical  channels  since  the  energy  is  averaged  over  the  entire 
objective  aperture.  The  beam  splitter  and  mean  intensity  de- 
tector will  be  added  if  it  proves  necessary. 

A block  diagram  of  the  receiver  electronics  is  shown 
in  pigure  11-5.  The  operation  of  the  signal  processing  elec- 
tronics is  controlled  by  a central  controller  inside  the  receiver. 
The  operational  sequence  is  initiated  by  external  stimulus  but 
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takes  place  under  internal  control. 

The  received  signal  is  detected  by  the  t IN  photo- 
diodes and  converted  from  current  to  voltage  by  the  transimped- 
ance amplifiers.  The  analog  processing  circuitry  for  each 
channel  consists  of  a variable  gain  amplifier,  an  integrator 
end  a sample  and  hold  unit.  DC  voltages  from  the  sample  and 
hold  are  converted  to  digital  signals  by  a 16  channel  multi- 
plexed 12  bit  A to  D converter.  This  digital  information  is 
stored  in  shift  registers  that  can  be  accessed  by  the  computer. 

An  operation  cycle  is  initiated  when  the  computer 
sends  out  a control  signal  requesting  DC  offset  data  from  the 
receiver.  The  integrators  are  gated  open  for  a precise  amount 
of  time  and  the  resulting  voltage  is  converted  and  stored. 

The  computer  then  reads  this  information  before  firing  the 
laser. 

After  the  laser  is  fired,  the  outgoing  pulse  is 
detected  by  a photodiode  external  to  the  receiver.  A signal 
from  this  detector  initiates  a timing  sequence  in  the  control 
electronics,  and  3tarts  the  counters  in  the  ranger.  After  a 
precise  interval  based  on  the  previous  range  reading  the  in- 
tegrators are  gated  open  to  capture  the  returning  pulse  and  the 
data  is  gated  and  stored.  The  returning  pulse  is  detected  by 
the  ranger  detector,  and  completes  the  time  of  flight  measure- 
ment that  indicates  range. 

One  millisecond  (this  time  delay  is  adjustable)  after 
the  first  'I’jlse,  the  second  pulse  is  transmitted,  and  the  process 
repeats  itself,  except  that  the  ranger  does  not  opente. 

After  the  information  from  both  pulses  is  stored,  it  is 
transferred  to  the  computer  where  the  data  is  processed  to  ob- 
t<j.l  variance,  covariance,  time  delayed  covariance  and  cross- 

wind. 
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The  signal  to  noise  ratio  for  the  receiver  was 
evaluated  using  the  techniques  described  in  an  earlier  report. 
Our  calculation  yielded  a signal  to  noise  ratio  of  31  db  with  a 
0.1  reflectivity  target  at  3 km  range  and  a receiver  gate  time 
of  3 microseconds.  The  lower  laser  energy  (100  millijoules) 
and  smaller  detector  aperture  in  the  1.06  system  were  partially 
offset  by  using  better  amplifiers  (wider  bandwidth  and  lower 
voltage  noise)  and  by  operating  the  photodiodes  fully  depleted 
for  minimum  capacitance. 

The  ranger  has  a resolution  of  2 meters  and  a calcu- 
lated root  mean  square  error  of  i 2 meters  with  an  0.1  re- 
flectivity target  at  3 km.  The  receiver  is .shown  in  Figures 
II-6  through  II-9. 

Two  modes  of  operation  are  designed  with  the  control 
software.  When  the  system  is  powered  up,  the  operator  selects 
the  number  of  double  pulse  cycles  (limited  to  less  than  41) 
that  are  to  be  used  for  one  sample.  This  preselected  number 
of  cycles  is  used  regardless  of  the  mode  of  operation  selected. 
In  the  normal  mode,  the  system  is  automatically  cycled  the  pre- 
set number  of  times  and  the  measured  wind  and  range  are  dis- 
played. Then  at  the  operator's  option  a hard  copy  of  the 
results  and/or  a paper  tape  of  the  results  or  raw  data  may  be 
obtained.  In  addition,  a video  display  of  the  raw  data  is 
provided  to  aid  the  operator  in  aligning  the  system  and  ver- 
ifying proper  operation. 

In  the  cw  mode  of  operation,  at  the  end  of  the  preset 
number  of  cycles,  the  results  are  punched  on  paper  tape  and 
then  the  system  is  automatically  cycled  again.  This  continues 
until  the  operator  terminates  the  cw  mode. 

The  gain  of  the  system  must  be  adjusted  to  take  into 
account  the  decrease  in  receive  i signal  strength  with  increasing 
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distance  to  the  target  and  differences  in  target  reflectivity. 
In  addition,  since  the  received  signal  is  a random  variable, 
consideration  must  be  given  to  the  expected  fluctuation  in 
signal  about  the  mean  value  so  as  not  to  limit  the  signal 
unnecessarily  and  distort  the  apparent  statistics. 

Our  previous  analysis^  assumes  that  the  received 
fields  are  Gaussian  distributed.  This  in  turn  leads  to  an 
exponential  probability  density  function  for  the  intensity 
which  can  be  written  as 


Pj(o) 


a 

<1> 


<l> 


Now  if  the  value  of  intensity  for  which  saturation  of  the  elec- 
tronics occurs  is  given  by  b,  then  the  fraction  of  samples  T 
on  the  average  for  which  the  receiver  will  be  saturated  is  given 
by 

b b 

T " J Pj(a)do  » 1 - e * 
o 

A reasonable  value  for  T is  0.95.  This  corresponds  to  b/^P  = 3. 
The  gain  can  then  be  set  such  that  the  mean  is  approximately  equal 
to  oue-third  of  the  peak  allowable  output.  This  same  principle 
can  be  used  to  set  the  gain  even  if  the  distribution  is  not 
exponential. 

The  system  gain  is  set  by  the  operator  through  the  use 
of  switches  on  the  control  panel.  Once  the  gain  has  been  properly 
set,  both  the  average  received  intensity  and  the  target  range  are 
monitored  by  the  computer  and  if  a significant  change  is  de- 
tected, the  operator  is  alerted.  It  is  anticipated  that  future 
versions  of  the  software  will  incorporate  automatic  gain  setting 
by  the  processor. 
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III.  ANALYSIS 

The  work  that  follows  is  based  on  the  extended  Huygens- 
Fresnel  formulation  and  includes  the  effects  of  the  turbulent 
atmosphere  on  the  laser  beam  as  it  propagates  to  the  target  and 
on  the  speckle  as  it  propagates  back  to  the  receiver.  Formu- 
lations are  developed  for  the  variance,  covariance,  time-delayed 
covariance,  slope  of  the  time-delayed  covariance  and  the  temporal 
power  spectral  density  of  the  received  intensity.  These  formu- 
lations include  the  effects  of  the  log-amplitude  covariance  as 
well  as  the  wave  structure  fui^tion  and  represent  a significant 
extension  of  previous  work. 

The  source,  target,  and  receiver  configuration  is  shown 
in  Figure  III— 1 . The  present  analysis  is  confined  to  the  case 
of  a TEMqq  laser  Illuminator.  It  is  assumed  that  the  source  and 
target  are  much  smaller  than  the  path  length  L,  that  the  distance 
between  the  receiver  and  source  is  much  smaller  than  the  path 
length,  and  that  the  outgoing  and  returning  radiation  experience 
independent  turbulence  regions.  The  last  assumption  is  thought 
to  be  correct  for  the  case  of  a perfectly  diffuse  target  when 
the  received  fields  are  jointly  Gaussian  (low  turbulence)  or 
when  phase  perturbation  is  the  dominant  turbulence  effect.  It 
may  not  be  a valid  assumption  when  the  log-amplitude  perturbation 
is  dominant  and  the  validity  of  the  last  assumption  under  this 
condition  remains  to  be  investigated. 

The  source  amplitude  distribution  can  be  written  as 


U (r)  ■ U exp 
o o 


/ _ri_  ikrA 

V V ' 2F  ) 


(Ill-l) 


where  a and  F are  the  characteristic  beam  radius  and  focal 
o 

length  respectively.  The  field  at  the  target  before  scattering 

from  the  target  then  is  written  from  the  extended  Huygens-Fresnel 
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principle  as 


2! 


Turbulence 


( I II— 2 ) 


where  describes  the  effects  of  the  random  medium  on  the  propa- 
gation of  a spherical  wave  from  the  source  to  the  target. 


The  field  at  the  receiver  is  written  by  reapplying  the 
Huygens-Fresnel  principle  to  the  field  at  the  target  and 


(III-3) 

where  U(p)  is  the  field  at  the  target  after  scattering  from  the 
target.  Utilizing  ( I II— 2)  and  (II1-3),  the  statistics  of  the 
received  intensity  can  be  formulated. 

Since  all  of  the  desired  statistics  can  be  derived 
from  the  time-delayed  correlation  function,  it  will  be  developed 
first.  The  time-delayed  correlation  function  for  the  intensity 
can  be  expressed  as 


B1(Pi.P2»t)  = <U(pi,o)U*(pi,o)U(p2,T)U*(p2,T)> 

4 

' (*)  ////  dp  ldP2dp3dpt4  <U(p  j ,o)U*(p2,o) 

U(p  3,T ) U*  (p  4 )p  expjj^  (Pi2  ~ P22  + P32  “ P42 
- 2 pi  * (p  1 - p2)  - 2 p2  ‘ (P 3 " P4))] 

h(p 1 , p2,  P3,  P4;  Pi»  p2;  t)  (111-4) 


where 
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H = <exp(i}>(pi , pj,  o)  + ij'*(Pi * P2,  o) 

+ ^(p2 » P 3 » t)  + i|/*(p2,  P4,  l)]> 

c exp  [-  j (D 12  - Dn  + Dl4  + D23  - D24  + D34) 

+ 2 C + 2 C 1 (III-5) 

X X J 

13  24 

and  where  the  wave  structure  function  and  the  log-amplitude 
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covariance  function  a’-e  given  by  ' 


2 Du  “ 2 <[*(v  pi' s)  - *(pi'  pr  tJ)]2> 
•hr^  S c„2  - 'F(tj  - ct) 

o 

r / 

+ (1  - t)  (pj  - pt)|  ' dt 


1 °°  / 

C * 0.132  tt2  k2  L f dt  C 2(t)  f dn  u 8 3 
Xij  J " J 

oo 

Sia'  ^t-L-L]  J0  [p|t(?j  - ?,)  - v(tj  - t,) 

+ (pj  - pi)(1  - '>!j 


(IIT-o) 


(111-7) 


and  where  p2  = p i ; p3  = p2;  p4  = p2;  tj  * t2  * o; 

and  t3  = t4  = i.  The  dummy  variable  t represents  the  distance 
from  the  source  to  the  field  point  normalized  by  the  total  path 
length  L.  Since  we  are  considering  propagation  from  the  trans- 
mitter to  the  target  and  also  from  the  target  to  the  receiver, 
care  must  be  taken  to  distinguish  between  this  normalized  path 
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length  variable  for  the  two  directions  of  propagation  if  the 
formulations  are  to  be  valid  in  the  case  of  non-uniform  turbu- 
lence C^2  and  crosswind  velocity  V and  if  a wind  weighting  function 
is  to  be  developed.  Consequently,  we  will  use  t to  indicate  the 
normalized  distance  from  the  target  to  the  receiver  and  t will 
be  replaced  by  1 - t when  propagation  is  f ’ .ora  the  transmitter 
to  the  target.  It  should  be  noted  that  (iII-6)  is  good  to 
order  n2  and  that  (III-7)  is  good  only  to  order  np  In  addition 
there  are  also  some  limitations  to  (III-5)  which  are  discussed 
in  Appendix  A where  it  is  derived. 

After  scattering  from  the  target,  the  fields  are 
20  21 

jointly  Gaussian.  * Consequently, 

<U(pi,  o)  U*(p2,  o)  U(p3,  t)U*(p4,  t)> 

■ <U(P1,  o)U*(p2,  o)><U(p3,  t)U*(p4,  Tp 
+ <U(p1,  o)U*(p4 , t)><U(p3,  t)U*'p2,  o)>  ( I I I— 8) 

For  the  case  of  a perfectly  diffuse  target,  (III-8)  becomes17,22-27 
<U(Pl,  o)U*(p2,  o)U(p3,  t)U*(p4,  t ) > 

<KPl*  °) ><  I (p  3,  Tp  6 (p  i - p2)  iS(p  j - p4) 

2 

j<U(P4,  o)U*(p4,  i)><U(p2.  t)U*(p2,  o)> 

<5(P  i - P„)  <S(p  3 - p?)  (III-9) 

Using  < I I I— 9)  in  (III-4)  and  performing  the  dp j and  dp3  inte- 
grations, the  time  delayed  correlation  function  becomes 

B (p,  r)  = — — lid  p2dp4  < I (p  o)>''l(P4  • 1 )>  H j(p?,  p4,  p ] , p?) 
u2l4 
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+ ~ ~ / f dp2dp4 

,k  - - 

ir  P • p 

e H2(P2 » 


<U(p4,  o)U*(p4 , t)><U(p2, 


Pt| , PI,  P2> 


i)U*(p2,o)> 


(III-10) 


where 

Hj  - H 

pi  * p2;  p3  ■ P4;  P2  “ pi*  P3  * P2 ; pn  ■*  P2 

4C  (p,  p,  t) 

- e x (III-ll) 


H2  - H 

pi  * p4;  p3  * P2 i P2  * pi;  p 3 ■ P2 ; p»>  ■ P2 


* exp  [20^(0,  p,  o)  - D^(p,  - p,  t) 

+ 2D^(p,  o,  t)  - D^(p,  p,  t)J  + 2C^(p,  " P»  t) 

+ 2Cx(p,  p,  i)J  (II 1—12 ) 

P * P2  " Pi 


p = p - p? 

4 

and  where  in  (II 1—6)  and  (II I—  7 ) » - p^  = p;  p^  - p^  = p;  and 

t - * t.  The  time  delayed  covariance  is  given  by 

Cj(p,  t)  = Bj(p,  t)  - <I(p>'2  (III— 13) 

where  from  previous  work^ 
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<kp)>  , _l  rdp  <kp)> 
nL2  J 


( 1 11—14 ) 


Utilizing  (III-13)  and  (III— 14)  in  (III-10)  the  time  delayed 
covariance  becomes 


4C  (p, 

CT(p,  t)  - CT  + C, — ' 1 J~  J-  X 


12  it2L4 

I(p2)><I(P4>> 


i rr  - - ( 4cx(p.  \ 
- JJ  dp2dp,  (e  - l) 


— tt 

v2LkJJ 


.k  - - 

*1  p • P 


dp2dp4  H2(p,  P,  t)  e 


<U(p4,  o)U*(p4,  t)>cU(p2»  t)U*(p2,  o)> 


Again  utilizing  the  extended  Huygens  Fresnel  principle 


(111-15) 


<U(p4,  o)U*(p4,  x)>  = U°  JJ  drj  dr2  exp 


rrl2  + r22' 


2a 


+ 1 2L  {l  ' f)  (r‘2  " r2?)  " 1 i p|*  ' (r!  " r?) 
*>exp[<J»(p4,  r!,  o)  + '(<*(p4 , r2 , t)]> 


(111—16) 


where 

<exp[sJ>(P4.  r j,  o)  + **(p4,  r?,  \)> 

= exp |~  | D^(o , r2  - r j , i )|  (II 1-17) 

Now  making  the  change  of  variable 
ri  - r?  = r 

r j + r>  = 2R 


2 h 


in  (III-16)  it  becomes 


<U(p4,  o)U*(p4,  t)> 


(A)  "o  !//dtd" exp  -(~ 


4R2 


ik 


(l  - £)  (r  • R)  - i£  e4  ' r - ± D^(o,  - r,  t) 
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(III-18) 


Performing  the  0 and  the  R integrations,  (III— 18)  becomes 
K 


<U(P4,  o)U*(P„,  t)>  - i Uo2  -f-  J dr 


exp 


4a  ' 
L.  o 


^ p4  • r - \ y°»  - r» 


[!f  (*-jOj 


(III-19) 

Using  (III— 19)  in  the  second  term  of  (HI— 15) , and  making  the 
change  of  variable 

P4  - P2  = p 

Pt,  + P 2 = 2R 
it  becomes 


CT  (P,  t) 


fk2V  2a  2\  ff 

/ JJ  dfl 


dr2  exp 


fr  )2  + rp2' 


4a 


\ V°»  ” r»*  T)  _ I V°*  " c2>  ~ T> 


-|f  (*-?)]  (r'2  + r^) 

(do  exp  ji£  p • 


P ” i2L  P ’ ^r2  ~ rl> 


(j  H2(p,  p, 
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J RdR,/*  d0R  exp  ["  *L  R ’ (r2  + rj)] 


(III— 21) 


Now  performing  the  6 and  the  rj  integrations,  (III-21)  becomes 

R 


CT  (p,  T) 


fe) 


<I><!  I dr2  exp  - 


2 


D , (o  , - r2 , i) 


-2r22[r  f (l  - ?)]  K 


(r2  + p) 


H , (p,  p,  i ) 


(III— 22) 


The  first  term  in  (III— 15)  will  now  be  evaluated.  The 
average  intensity  can  be  obtained  by  letting  t ■ o in  (III— 19) 
a, id  then  performing  the  6r  integration.  It  becomes 

,2  U 2a  2 f . 2 

■ k)  -V-  / rdr  jo(l  p-r)  exp  • 7T 

L o 

- i D^io,  - o)  - [i  f (l  - i)]  r;’  (111-23) 


Using  (II 1—23 ) in  the  first  term  of  (III-15),  and  making  the 
change  of  variable 

p2  - Pi,  = P 

Pi*  + P2  = 2R 


it  becomes 
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ri,  o) 


1 D*<°*  - r2- 0>  - (r  f l1  - f))  ('*' + r^) 


(II 1—24 ) 


Now  performing  the  0 , R,  and  r integrations,  (III  — 24 ) becomes 

K 

V*’  ’>  ‘ if  (l)ff  r2dr*  dd  Jo(r  r2p)(e  x<l>’ C’  l) 


exp 


2£_ 

2a  2 
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- Df(o.  r?,  o)  - j(£f  (l  - i)) 


(II 1—25) 


This  completes  the  derivation  of  the  time  delayed  covariance 
function.  It  is  represented  by  the  sum  of  equations  (II 1—22) 
and  (111-25). 

There  are  four  statistical  quantities  that  are  of  in- 
terest and  can  be  derived  from  the  time  delayed  covariance 
function:  the  variance;  covariance;  power  spectral  density;  and 
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the  slope  of  the  time  delayed  covariance  function  at  zero  time 
delay.  The  last  quantity  is  of  special  interest  since  it  is  pro- 
portional to  the  path  averaged  crosswind.  The  other  statistical 
quantities  are  needed  for  overall  system  design. 

The  variance  can  be  obtained  from  the  time  delayed  co- 
variance  by  setting  p * o and  r = o.  ‘Under  these  conditions, 
the  0p  integration  can  be  performed  and  the  variance  is  given  by 


<I>‘ 


- r2‘ 


-D^(o,  r2,  o) 


(i:)//r2dr2pdp  exr 

i i (*  - D)]  p<0>- 1]  »>) 


2a  2 
• o 


(111-26) 


The  covariance  can  be  obtained  from  the  time  delayed 

covariance  by  setting  t ■ o.  This  allows  the  0 integration  to 

r2 

be  performed  in  the  second  term  and 

CT(p)  - C.  (p,  o)  H-  CT  (p,  o) 
i 1 1 i2 


where 


Cj2(P*  o) 


Hr  (riff 


r2dr2  exp 


, 2a 

i-  o 


2r,2 
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[l  f (>  - 1) 

mm 

f - ‘ 

I dp  e 


2 <l> 
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D (o,  - r2,  o) 


- *1  P ' p 


Jo  (l  pr?)  H2(p’  P»  °> 


(II 1-27) 
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and  Cj^(p,  o)  is  given  by  (III-25)  with  t = o. 

The  temporal  power  spectral  density  is  the  Fourier 
transform  of  the  time  delayed  covariance  function  evaluated  at 
p * o.  Since  no  additional  integrals  can  be  evaluated  it  will 
not  be  written  in  detail,  but  it  can  easily  be  obtained  from 

OD 

sn(w)  = y*  e ^WT  jCj  (o,  t)  + C^(o,  T)]dT  (111-28) 

_ oo 

and  equations  (II 1-22 ) and  ( I 11-25) . 

The  slope  of  the  time  delayed  covariance  function  is 
given  by 

MI(p)  " I?  Cl(p*  T) 

T = O 


57  ci,(p'  T), 


+ 57  ci2(p’  T) 


(1 1. 1-2  9) 


where  the  two  derivatives  with  respect  to  t can  be  obtained  from 
( I 11—22)  and  (III-25).  Using  ( I 11-25) , 


57  ci  !<p'  ’> 


(!/4//r 2dr^p  Jo(c 


4C  (p,  p , o) 
„ X 


T = O 


- D (o.  r2»  °)  a 

6 f (r2 } fT‘  Cx(p.  P,  t) 


T = O 


(111-30) 
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where 


S7  cx(p*  P*  T) 


f dt  c 2 V • c (p,  p) 

* » X 


(III— 31) 


T ■ O 


C (p,  p)  - 0.132  IT2  k2  L 


pt  + (1  - t)pj 


/ du  u': 


Sln2(~--^^— t-)-)  J ,(ujpt  + (1  - t)p|)  ( 1 11—32) 


f (r, ) * exp 


-2(rf  v1-^))^' 


From  ( I 11—22) . the  second  term  in  (111-29;  becomes 


vF-  ,>|-^(t)7‘;»  f<r2)/d?e 


k~  — — 

i-p  • ip  + r?) 


Hr,  (p , p,  t) 


■ t(r?)/ dire ' 


i-p  • (p  + r2) 


-D,(o,  r?,  o) 
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H?(p,  P,  o) J dt  C^2 


•['V  <°»  - r2)  + \ V (P.  -P)  - % <P»  °> 

+ £ V <p,  p)  + 2 cx  (p,  - P)  + 2Cx  (?,  p)j  ( III— 33) 

where  use  has  been  made  of 


r,' D* (p' t,>T) 
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dt  C„J  v • v<p- »> 


t * o 


and 

D <?,  P>  * -9*70  Lk"  •?£' 1/3 

V |pt  + (1  - t)p  I 


(III-34) 


(II 1—35) 


The  0r2  integration  can  now  be  performed  on  all  but 
one  term  and  the  slope  becomes 
2 . 1 

(IT  1—36) 


v?>  - /k 
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dt  C 2 V • U 
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where  the  vector  wind  weighting  function  W is  given  by 

W =jfy*  r2^r2  dp  f(r2) 


-D  (o,  r2,  o) 


- — — *” 

- 4C  (p,  p,  o)_  _ _ _ _ if  P * P / . 

4e  x C^Cp,  (>)  + H2(p’  P’  °^e  \2 

- D (p,  o)  + £ 0 (p,  p)  + 2C  (p,  ~ p)  + 2C  (p,  p>) 

X.  X_  / 


33 


f (r2)dp 


(p  + r2)  - D (o,  r2,  o) 


H2(p,  p,  o)  (o,  - r2) 


(III-37) 


The  results  of  this  section  are  summarized  in  Table 
III-l  which  lists  the  pertinent  equation  numbers  for  such  sta- 
tistical quantity  that  has  been  developed. 


Table  III-l.  Equation  Summary 


Statistic 

Time  Delayed  Covariance 

Variance 

Covariance 

Power  Spectral  Density 
Slope 


Pertinent  Equations 

13,  22,  23,  6,  7,  12 

26,  6,  7 

27,  25,  6,  7.  12 

28,  22,  25,  6,  7.  12 

36.  37.  6,  7,  12,  32,  35 


Numerical  evaluation  of  these  equations  Is  considered  in 
Section  IV, 

A physical  interpretation  of  the  variance  and  covariance 

reveals  the  mechanisms  Involved  in  speckle  propagation  through 

turbulence  in  a particularly  illuminating  manner. 

The  variance,  (111-26),  heuristicaily  may  be  expressed  as 

6Y<>  ? 

X 

o ? - 2 »•  - 1 (I  11-38) 

n 

which  in  cases  where  4 i u ■'  • • l can  be  written 

A 

o 2 * 1 + 8yo  7 (11 1-39) 

I X 

n 

In  these  expressions,  <’  ■ is  the  log  amplitude  variance  for  a 
point  source  propagating  over  the  target-receiver  path,  and 
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Y (<  1)  is  the  "source  smoothing  factor"  associated  with  the 
nonzero  spot  size  (S)  on  the  target.  Mathematically,  y relates 


to  the  integration  over  (p)  in  (III-26),  with  the  spot  size 

(S  = or  coupled  through  r.  For  a small  transmitter 

ka0  kp0 

(or  defocused  spot),  or  strong  turbulence,  the  target  spot  will 

be  S > /i/k  and  y will  be  less  than  unity.  In  the  limit  of  strong 

turbulence,  y =»  0 while  for  weak  turbulence,  o '*  = 0;  either 

extreme  results  in  a unity  variance. 

A quantitative  estimate  of  y may  be  obtained  as  follows. 
28 

From  published  results  for  the  variance  with  an  extended  in- 
coherent source,  we  know  that  in  weak  scattering 
, o v-7/3 


(/dk) 


( 1 11-40) 


The  smoothing  in  the  strong  scattering  case  is  conjecre**«d  to  be 


(111-41  ) 


This  source  or  target-spot  smoothing  is  analogous  to  finite- 
receiver  smoothing. 

The  expressions  corresponding  to  (111-38)  and  ( 1 11—39) 
for  an  incoherent  source  are 


4yu  ^ 
X 


°1  ‘ = ° 
n 


(111-42) 


o ^ 4yo 

!n  X 


(I  1 1-4  1) 


It  may  be  seen  that  the  unity  baseline  and  half  the  remaining 
term  in  (111-39)  arise  t rum  the  source  coherence  or  speckle 
phenomenon,  whether  controlled  by  the  target  (weak  turbulence) 
or  atmosphere  (strong  turbulence).  Alternatively,  we  may  observe 
that  the  unity  term  arises  strictly  1 rom  phase  randomization. 
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with  associated  exponential  irradiance  statistics,  while  the 
remaining  term  represents  the  effects  of  nonzero  amplitude  per- 
turbations, with  associated  log  normal  statistics. 

To  further  illustrate  these  considerations,  let  us 

consider  the  large,  focused  transmitter  condition  (a  >>  /L/k) • 

o 

We  first  consider  the  weak  turbulence  condition  (p  > a > / L/k , 

o o 

in  which  the  transmitter  is  able  to  operate  in  an  essentially 

diffraction-limited  manner.  The  target-spot  is  small,  i.e. 

such  that  the  receiver  is  in  the  far-field  of  the  spot,  and 

the  situation  is  much  like  that  of  a point  source  (a  2 l,y  = 1). 

X 

There  is  however  one  major  difference  from  the  point-source  case: 
we  assume  a good  statistical  (ensemble)  average  of  the  field, 
which  as  a practical  matter  means  that  the  target  Is  moved  slowly 
so  that  the  (large)  speckles  sweep  past  the  receiver.  (If  the 
transmitter,  target,  and  receiver  were  rigidly  fixed,  the  re- 
ceiver would  see  part  of  one  speckle,  and  the  turbulence  fluctu- 
ations would  be  that  of  a point  source,  with  purely  log  normal 

v2  o 

statistics  and  a variance  of  e A - ! % 4o^)»  The  resultant 
variance  Is  essentially  unity,  with  small  turbulence-caused 
fluctuations  superimposed  on  the  basic  speckles  that  would  be 
present  in  the  absence  of  the  atmosphere.  The  irradiance  sta- 
tistics will  be  nearly  exponential. 

We  now  let  the  turbulence  strength  grow  until  "saturation" 

29 

or  multiple  scattering  sets  in.  The  peak  value  of  o i is  known 
to  be  approximately  0.6,  occurring  when  p^  is  on  the  order  of  /i/k. 

The  target  spot  will  no  longer  be  diffraction  limited  (a  p 

o o 

/l./k),  and  hence  y will  be  less  than  unity.  The  exact  value 
of  the  peak  variance  can  only  be  determined  i rom  a lull  numerical 
evaluation  of  (111—26) , including  the  detailed  effects  of  satu- 
ration on  C (p);  the  statistics  will  be  a composite  of  log  normal 

X 

and  exponential  distributions.  As  the  turbulence  strength  grows 
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further,  a 2 decreases  ("super-saturation")  to  the  limit  5:5  0.17; 

X 

also  the  target-spot  grows  indefinitely,  so  that  y decreases  to 
zero  *n  the  limit.  The  variance  therefore  decreases  back  to 
unity,  with  exponential  statistics. 

We  now  consider  the  covariance  case,  again  for  a 
focused  transmitter.  For  the  coherent  term  given  by  (II 1—27 ) 

if  C is  neglected  and  either  a or  p is  much  smaller  than  the 

X o o 

other  and  also  smaller  than  /L/k,  the  corresponding  limitation  on 
r in  the  integrand  couples  through  the  Bessel  term  and  exponential 
to  imply  a similar  limitation  on  p.  This  gives  a closed-form 
covariance  result  which  is  that  obtained  with  the  assumption  of 
joint  gaussian  field  statistics^’ ^ and 

~p2/2a  2 - 4(p/p  )5/3 

CT  (p)  - e 0 0 (111-44) 

n 

Hence  the  covariance  arising  from  the  coherent  term  is  generally 

exponential^  with  a scale  equal  to  the  smaller  of  or  pQ 

(target  and  atmospheric  speckles  respectively).  In  the  event 

that  a and  p are  comparable,  the  complete  expression  must  be 
o o 

used  and  the  field  statistics  are  marginally  (but  not  joint) 
gaussian. 16,17  in  the  event  that  C is  non-negiigible,  the  co- 

A 

variance  scale  relating  to  the  coherent  term  is  largely  unaffected. 

We  now  consider  the  incoherent  term  (111-25).  If  the  target 

spot  is  small  compared  to  /L/k,  i.e.  u and  pq  /L/k,  the  large 

range  of  r in  the  integrand  will  couple  through  Jq  to  hold  p small 

in  C . Hence  the  covariance  scale  will  be  determined  by  C (p,0) 

X X 

and  will  be  /L/k  as  expected  for  a small  target  spot.  However, 

for  a large  target  spot  (S),  the  p range  in  the  integrand  will 

be  large  and  the  p-scale  in  C will  be  "smeared"  to  the  order  of  S. 

This  effect  is  consistent  with  the  incoherent-source  results  of 


Fef.  8,  and  can  be  argued  from  simple  geometry:  for  two  receivers 

separated  by  a small  distance  relative  to  the  source  size,  the 
source-turbulence  configuration  over  the  paths  to  each  respective 
receiver  will  be  highly  correlated.  The  statistics  associated 
with  the  incoherent  term  are  log  normal,  although  the  log  amplitude 
is  not  jointly  normal  except  for  weak  fluctuations. 
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IV.  NUMERICAL  EVALUATION 

The  theoretical  formulations  for  the  statistics  of  the 
received  intensity  must  be  reduced  to  numbers  as  required  for 
system  evaluation  and  design.  Standard  numerical  integration 
techniques  proved  too  time  consuming  due  to  the  three-fold 
integrations  required,  and  Monte-Carlo  techniques  would  not  work 
due  to  the  form  of  the  integrand  and  the  interdependency  of 
integration  limits.  However,  a major  analytical  breakthrough 
has  been  achieved  which  greatly  simplifies  the  formulations 
while  retaining  a complete  description  of  the  mechanisms  in- 
volved. We  are  now  able  to  reduce  the  formulations  to  numbers 
for  comparison  with  the  cw  experimental  results  and  for  design 
purposes. 

The  variance  and  covariance  have  the  following  form 
for  the  case  of  uniform  turbulence 


*cx(p) 


1>2  (l)2  //  rdr  0d0  [2e  X - 1]Jo(!pr) 


f?(r)  (IV-1) 


0 o 


C (p)  = C (p)  + C (p) 
I 1 1 12 


( IV— 2 ) 
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(IV— 4 ) 


drdp  rp  f j (p , p)  rp)  f2(r) 

5/3 


f i(p , P)  = exp 
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r i 


/ dt|pt  + (1  - t)p|  ' + / dt | pt  - (1  - t)p | 


_ 5/3 


L o 


+ 2 C (p,  p)  + 2 C (p,  - p) 
X X 


f 2(r)  - exp 


L ' o 


(IV-5) 


( I V— 6 ) 


Each  of  the  integrands  contains  the  term  rJ^r  rp^  f2(r).  0ur 
approach  is  to  expand  f2(r)  in  a Fourier-Bessel  series  and  then 

make  use  of 

00 

26 (a  - B) 


/26(a  - 

rJo(ar)J0(6r)dr  = T7T 

o 

to  simplify  the  integration.  The  series  is  given  by 

»*<*>  J„  & r) 


( IV-7 ) 


(IV— 8) 


where 


m a 2 ,2 


A^J  (P  ) 
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/ *f«*>  Jo  (f  *) 


dx 


( 1V-9) 


and 
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J (pJ  = 0 

o in 


This  expansion  can  be  done  because  f2(r)  becomes  negligible  for 
some  value  of  r that  is  easily  calculated  and  does  not  depend  on 
the  other  integration  variables.  In  addition,  the  shape  of  f2(r) 
allows  it  to  be  represented  by  only  a few  terms  in  the  series. 

Our  experience  indicates  that  for  problems  of  interest  to  this 
program,  around  six  terms  in  the  series  yields  accurate  results. 

Using  (IV-7)  - (IV-9),  the  variance  and  covariance 


become 


<l>2  £ b 


m 


ra 


1 


(IV-10) 


(IV-11) 


where  0 is  the  angle  between  P and  p.  The  net  effect  of  the 
above  work  is  in  the  case  of  the  variance  to  reduce  a twofold 
integration  to  a set  of  well  defined  onefold  Integrations; 
and  in  the  case  of  the  covariance  to  reduce  a rather  involved 
threefold  integration  with  interacting  limits  to  a well  defined 
onefold  integration  plus  a set  of  well  defined  onefold  inte- 
gral ions. 

Equations  (IV-10)  and  ( IV— 11)  ich  contain  the  log- 
amplitude  covariance  function  given  to  >rder  nj  by^ 
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CX(P,  p)  - 0.  132  7T2k2Cn2  J dt  j du  8/3  sln'?(“2't~  ~ -t)L) 

0 O 

JQ(u|pt  + p(l  - t)|)  ( I V— 12) 

Evaluation  of  ( IV— 12 ) Is  now  the  limiting  factor  in  generating 

numbers  from  the  variance  and  covariance  formulations.  For  the 

variance,  p =*  0 and  (IV-12)  can  be  represented  by  a six  term 

series.  This  has  allowed  us  to  develop  an  efficient  program 

that  readily  computes  the  variance.  This  program  was  used  to 

generate  the  analytic,  unsaturated  variance  curves  shown  in 

Section  V which  discusses  the  cw  experimental  results. 

For  the  covariance,  p is  not  zero  and  (IV-12)  is  not 

30 

easily  reduced  to  a simple  form.  The  technique  used  by  Fried 
for  the  case  p =>  0 unfortunately  does  not  work  when  p j 0 and  for 
the  parameter  values  of  interest.  Several  other  expansion  tech- 
niques have  been  attempted,  but  so  far  we  have  not  been  able  to 
reduce  (IV-12)  to  a more  tenable  form.  We  are  still  working  to 
resolve  this  problem.  Meanwhile,  a twofold  Gaussian  quadrature 
numerical  integration  is  being  used  to  evaluate  (IV-12).  It  is 
relatively  time  consuming  because  for  problems  of  practical 
interest,  many  cycles  of  the  SIN  squared  term  must  be  included  in 
the  range  of  integration.  However,  we  have  been  able  to  generate 
some  numerical  values  for  the  unsaturated  covariance  as  shown  in 
Section  V. 

Equation  (IV-12)  holds  only  for  non-saturated  turbu- 
lence conditions.  For  saturated  turbulence  there  is  presently 
no  two  source  log-amplitude  covariance  formulation  available 

for  use  in  ( I V— 11)  and  one  must  be  developed.  However,  Clifford 
31 

and  Yura  have  developed  a single  source  (p  = 0)  formulation. 
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For  strongly  saturated  turbulence  it  is  given  by 

] 00 

Cx(p)  « 0.42  a 5/2 J du  J dq  e"q  q'3/5  J (IV-13) 

oo 

where  the  parameter  a is  determined  from  a knowledge  of  the 
saturated  value  of  the  normalized  variance  of  the  intensity  for 
a spherical  wave.  Using  unity  as  the  saturated  limit,  we  have 

C (0)  - .17329  - 0.42  a5/3  / du  / dqe"q 
X *o  Jo 

and  the  value  for  a is  * 0.588.  Equation  (IV-13)  then  becomes 

1 

du  J /s.693  uq  3/5  2-\  (IV-14) 

°\  V 

Equation  (IV-14)  has  been  evaluated  using  standard  numerical 
integration  techniques  and  used  in  conjunction  with  (IV-10)  to 
evaluate  the  variance  under  saturated  conditions.  The  results 
are  shown  in  Section  V. 

Equations  (IV-10)  and  ( IV— l 1 ) should  be  valid  for  all 
turbulence  conditions.  However,  this  is  only  true  if  the  right 
form  for  the  log-amplitude  covariance  function  and  the  wave 
structure  function  are  used.  There  is  currently  no  simple  formu- 
lation available  for  these  quantities  valid  in  the  t.ansition 
region. 


Ve) 


17329 


/ d’e'<i  / 
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V. 


CW  EXPERIMENT 


The  experimental  observation  and  measurement  of  sta- 
tistical properties  of  the  intensity  after  reflection  from  a 
diffuse  or  quasi-dif fuse  target  and  propagation  through  atmos- 
pheric turbulence  are  presented  in  this  section.  The  purpose 
of  the  experiments,  definition  of  the  parameters  and  measure- 
ments of  interest  and  the  experimental  design  will  be  discussed. 
In  addition,  the  results  of  experiments  performed  to  date,  in- 
cluding the  relationship  to  recent  analytical  and  numerical 
work  will  be  r.esented. 

The  observation  of  scattering  from  a coherently 
illuminated  diffuse  target,  is  simply  illustrated  in  Figure  III-l 
showing  the  source  (laser),  target  and  the  receiver  coplanar  with 
the  source.  The  source  and  receiver  are  separated  by  a distance 
that  is  large  compared  to  the  scales  of  interest  in  the  received 
intensity  such  as  aQ  (transmitter  radius),  /l,/k  , and  pQ,  and 
very  small  compared  to  the  source-receiver  and  target  separation. 
The  experiments  described  below  Involve  the  measurement  of  sta- 
tistical quantities  discussed  analytically  in  previous  Sections 
and  other  reports.  These  include  the  probability  distribution 
function  of  the  irradiance,  he  normalized  variance,  and  the 
spatial  covariance  function.  In  situ  atmospheric  parameters 
including  temperature,  pressure,  wind  velocity  and  C^2  were  also 
monitored  during  the  data  collection. 

The  experimental  system  is  schematically  illustrated 
in  Figure  V-l.  The  laser  is  an  argon  ion  laser  operating  at 
0.488  Urn  and  yielding  approximately  700  milliwatts  of  power. 

The  transmitted  beam  is  amplitude  modulated  at  100  kHz  using  a 
photo  elastic  modulator  and  a calcite  polarizer.  Scotchlite 
paper  is  used  as  the  target  material  because  of  its  high  gain  in 
comparison  with  a true  Lambertian  reflector.  Laboratory 
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experiments  were  conducted  verifying  the  'diffuse'  nature  of 
the  reflection  from  a scotchlite  target.  This  requires  that  the 
spot  size  on  the  target  be  large  compared  to  the  bead  size  of 
the  target  material  and  was  always  true  in  the  atmospheric 
experiments  conducted.  The  receivers  consist  of  two  S-20  photo- 
multiplier tubes  followed  by  10  kHz  bandpass  filters,  synchronous 
demodulators,  and  Bessel  low  pass  filters  with  switchable  band- 
widths  hetween  800  Hz  and  2.8  kHz.  All  of  the  statistical  data 
reduction  was  done  on  a digital  computer. 

Initial  experiments  using  a HeNe  laser  at  0.6328  microns 
over  a short  path  indicated  excellent  agreement  between  the  ex- 
perimental results  and  the  theory  under  low  turbulence  conditions. 

As  predicted  by  the  theory,  a normalized  variance  equal  to 
unity  was  observed.  Initial  experiments  with  the  argon  ion  laser 
however  produced  normalized  variance  values  of  the  order  of  0.3. 

This  was  significantly  lower  than  expected.  However,  the  co- 
variance  results  were  substantially  in  agreement  with  the  theory. 

The  low  normalized  variance  was  related  to  the  bandwidth  of  the 
argon  ion  laser  which  is  of  the  order  of  10  gigahertz.  This 
translates  to  a coherence  length  of  approximately  3 cm.  Because 
the  il’umin.ited  spot  size  at  the  target  is  of  the  order  of 
centimeters,  the  very  short  coherence  length  of  the  laser  will 
cause  a decrease  in  the  ability  of  light  from  different  points  on 
the  target  to  interfere  with  a resultant  decrease  in  the  observed 
normalized  variance.  Consequently,  the  laser  was  modified  to 
include  a temperature  stabilized  etalon  which  limits  the  band- 
width to  that  of  a single  mode.  The  results  for  the  normalized 
variance  later  in  this  section  indicate  that  this  modification 
does  solve  the  coherence  problem  observed  in  the  initial  experiments. 

A further  recent  modification  to  th“  experimental  system 
consisted  of  the  addition  of  a telescope  to  the  front  of  the 
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receivers.  The  actual  size  of  the  receivers  is  2 mm  and  the 
minimum  separation  is  A. 5 mm.  Under  high  turbulence  conditions, 

Pq  could  become  small  enough  that  the  covariance  function 
measurements  would  be  smoothed  and  increased  because  of  tne  size 
of  the  detectors,  and  the  variance  measurements  would  be  de- 
creased by  aperture  smoothing  over  several  speckles.  Two  tele- 
scopes were  constructed  which  could  be  set  onto  the  receiver  to 
reduce  the  effective  receiver  size  and  decrease  the  effective  mini- 
mum receiver  separation.  The  effects  of  the  telescopes  are 
summarized  in  Table  V-l. 

Receiver  Minimum  Relative 
Size  Separation  Received 


Magnification 

(mm) 

(mm) 

Power 

no  telescope 

1 

2.0 

A. 5 

1.0 

with  telescope 

1.7 

1.18 

2.65 

0.35 

with  telescope 

3. 1 

0.38 

1.A5 

0.10 

Table  V-l  Effect  of  Receiver  Telescopes 

The  use  of  the  telescopes  of  course  decreased  the  average  power 
received  and  thus  a compromise  was  required  between  the  obser- 
vation of  desired  conditions  and  the  limitations  of  the  receiver 
to  maintain  suitable  signal  to  noise  ratios. 

Each  experiment  consisted  of  a series  of  recordings  of 
the  irradiance  at  two  points  in  the  receiver  plane  separated  by 
a variable  distance.  From  these  recordings,  probaoiiity  distri- 
bution functions  (PDF)  were  computed  from  a single  channel,  and 
covariance  curves  were  constructed  by  measuring  the  correlation 
coefficient  between  the  two  signals  tor  different  separations. 

Figures  V-2  and  V-3  show  typical  PDF  calculations 
from  the  speckle  measurements  for  a focused  transmitter  and 


Measured  Probability  Distribution  of  Irrndiance, 
Weak  Turbulence  (l.  = 500  m,  o ' ’ - 1.02). 

*N 
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path  length  (L)  = 500  m.  The  first  of  these  figures  is  a low 

turbulence  condition  p (**  10  cm)  > a (1.3  cm)  > /h/k  (0.62  cm) 

o o 

and  thus  target  speckles  would  be  expected  to  predominate.  The 
expected  PDF  under  these  conditions  is  an  exponential  distri- 
bution.^ A corresponding  exponential  distribution  is  plotted 
in  the  figure  and  the  figure  is  well  within  the  scatter  of  the 
experimental  points.  The  calculated  normalized  variance  is  1.02 
as  compared  to  the  value  of  unity  predicted  by  the  theory. 

Figure  V-3  shows  the  results  for  a higher  turoulence  case  where 
the  log-amplitude  perturbation  is  important.  The  PDF  in  this 
case  shows  a marked  deviation  from  exponential  at  low  and  high 
values  of  received  intensity.  This  was  expected  since  the 
total  perturbation  in  this  case  is  the  result  of  the  combiner  ion 
of  two  processes  one  of  which  is  log-normally  distributed.  The 
measured  normalized  variance  in  this  case  is  1.2. 

Measurements  of  the  normalized  variance  o 2 over  a 

N 

wide  range  of  turbulence  levels  are  shown  in  V-4.  At  very 

low  turbulence  levels  (C^2  < !0“15  m~2^  *)  and  high  turbulence 

levels  (C^  2*  2 x 10-l!  m-?/!),  •’  is  generally  near  unity 

as  predicted  hy  both  tin*  simple  theory  and  the  phase  dominance 

theory.  The  intermediate  region  shows  an  increase  In  u ? above 

N 

unity  as  predicted  by  the  complete  theory.  AlKo  shown  in  this 
figure  are  the  results  of  numerical  evaluations  of  2 including 
terms  containing  C (p)  (log-amplitude  covariance  function).  The 

A 

left-hand  curve  was  calculated  using  a form  of  C (p)  valid  under 
low  turbulence  (unsaturated)  conditions.  The  right-hand  curve 

was  similarly  derived  using  a high  turbulence  model  for  C (p). 

X 

While  not  showing  particularly  good  quantitative  agreement  at 
intermediate  turbulence  levels,  the  theoretical  formulation  coupled 
with  the  experimental  results  have  provided  us  a good  understanding 
of  the  physical  mechanisms  involved  in  this  problem. 
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Normalized  Variance  for  a Focused  Transmitter  vs*  Log  Amplitude  Variance 
for  Point  Source.  (Etalon  Inserted)  L = 500  m. 


Results  of  typical  measurements  of  the  covariance 
functions  are  shown  in  Figures  V-5  through  V-10  for  a range  of 
turbulence  conditions  and  for  two  values  of  pathlength.  Low 
turbulence  conditions,  such  as  seen  in  Figures  V-5  through  V-8 
show  broad  covariance  curves  described  well  by  the  Joint 
Gaussian  theory.  As  the  turbulence  level  increases  and  p 
approaches  a^,  the  covariance  curves  narrow  as  seen  in  Figures 
V-9  and  V-10  and  differ  considerably  from  that  predicted  by  the 
joint  Gaussian  theory.  However,  as  shown  in  Figure  V-10,  there 
is  excellent  agreement  between  these  measurements  and  the 
complete  theory. 

To  date,  excellent  progress  has  been  made  in  verifying 
our  theoretical  formulations  with  cw  experimental  measurements. 
The  next  step  in  this  effort  is  to  verify  the  time  delayed 
statistics  from  which  the  crosswind  is  derived. 


e for  Focused  Transmitter,  Weak  turbulence.  (L  - 910 
cm,  / L/k  = 8. A ram,  a = 1.35  cm.) 


p(mm) 


Figure  V-10  Covariance  for  Focused  Transmitter,  Strong 

Turbulence.  (L  = 500  m,  p = 8. A mm,  /L/k  = 6.2  mm, 
a =1.35  cm.) 
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APPENDIX  A.  GENERALIZED  POINT  SOURCE  MUTUAL  COHERENCE  FUNCTION 


* 


The  extended  Huygens-Fresnel^  ^ approach  is  frequently 
utilized  for  determining  the  moments  of  an  electromagnetic  field 
after  propagation  through  a turbulent  medium.  Crucial  to  the 
use  of  this  method  is  the  availability  of  an  appropriate  point 
source  mutual  coherence  function.  It  is  the  purpose  of  this 
Appendix  to  develop  a generalized  n^1  order  mutual  coherence 
function  that  can  be  used  in  the  extended  Huygens-Fresnel  inte- 
gral to  generate  moments  of  any  order. 

The  generalized  nt*1  order  spherical  wave  mutual  coherence 
19 

function  is  defined  as 

H(p  1 1 P 2 * P2n;  i’l*  P?***P2n>  * ‘-exp[ij/(pi,  Pi)  + <!'*((>?.  P;) 

+ *(P?n-i,  P?n-l)  + '!'*<P2n,  P?n>)'  U> 

where  i|»  = x + i«l> * corresponds  to  a source  point  in  the  trans- 
mitter plane  and  p^  corresponds  to  an  observation  point  in  the 
leceiver  plane.  In  the  work  that  follows,  j>  will  be  treated  as 
a complex,  jointly  normal  random  variable.  With  ^ normally 

o 'X') 

distributed,  then  v x'  = and  * - 0.  It  should  be  noted 

that  this  is  not  an  additional  assumption,  but  the  result  of 
assuming  1*1  to  be  jointly  normal.  The  implication  of  tins 
assumption  on  <i>  is  discussed  in  the  summary  paragraph.  Conse- 
quently the  expected  value  in  (1)  can  be  broken  into  two  parts 
which  will  be  evaluated  separately. 

H - <exp{x(p i»  Pi)  + + X<P2n,  P2n)l>  <exp[i<KPi,  Pi) 

-14(P2»  P2)  + + i(KP2n-l.  P2n-i>  ' i<^^2n,  P2n))>  <2> 


6 4 


*3  - x(P3i  Pj) 


l<Xl> 


,<X2n> 


is  the  covariance  matrix  for  the  log-amplitude,  the  underbars 
denote  matrices  and  T indicates  the  transpose.  The  second  ex- 
pected value  in  (2)  is  given  by 
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and  is  the  covariance  matrix  for  the  phase  perturbation  term. 
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Equations  (3)  and  (4)  can  now  be  combined  to  form  the  mutual 
coherence  function.  However,  first  we  will  express  (3)  in  terms 
of  the  log-amplitude  structure  function  D and  covariance 


Ai1 

function  C and  will  express  (4)  in  termsJ  of  the  phase 
Xii 

structure  J function  D, 


where 


’xy  ' <IX<V  PJ)1!>'2(°X2  ■ \ j 

^ ^ ■ *<Pj.  Pj)32>-2(%2  • 


and 
C 


x “ <Ix(p1.  Pt)  - <X>]  tx(Pj»  Pj)  “ <X>J> 
Equation  (4)  can  then  be  rewritten  as 


- f"2n  2n  i+j  2n  2n  i+j 

•f  h I (-1)  O 2 + Z Z [-1]  fA  - o 2] 
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, 2n  2n-l  i+j 
+ f Z Z (-1)  D. 
j-i+1  i-l  *ij 


(5) 


and  equation  (3)  as 


-2n  o 2 
X 

e 


1 

2 


2n  2n 
Z Z 
J-l  i-l 


i+j 

(-1) 
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(6) 


2n  2n 


i+j 


+ f £ £ [1  + (-1)  1 A 

j-1  i-1  xij 


2n  2n-l  i+j  2n  2n-l  i+j 

1 £ £ (-1)  D + £ £ [1  + (-1)  ]C 


2 j-i+1  i-1 


lij  j-i+1  i-1 


lij 


Using  (5)  and  (6)  in  (2)  the  generalized  mutual  coherence 
function  becomes 


2n  2n-l  i+j 


2n  2n-l 


i+j 


£ (-1)  D,  + £ £ [1  + (-1)  ]C 


j-i+1  i-1 


'ij  j-i+1  i-1 


H - f! 


(7) 


where 


D.  * D + D. 

*u  xu  ♦« 


The  two  source  wave  structure  function  and  the  two  source  log- 
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amplitude  covariance  functions  are  given  by 

1 

D*iJ^i*  ***  Pj’  ^ " 2’91  C°2  Zk2/ 


+ (1 


5/3 


- t)(pj  - Pt)|  dt 


(8) 


U 

CX^Pi’  Pi’  Pj’  Pj^  " 0,132  "2k2  f dsCn2(8) 

o 

m 

■fju  » 8/3  .in2  ["■££'*>]  J0  [u|(Pj  - pt)  l 

+ <?j-P1)U-f)l] 


(9) 


6b 


It  should  be  noted  that  (8)  is  in  general  good  to  second 

order  in  ni  and  for  the  case  where  ip  is  jointly  normal  it  is 

19 

exact  to  all  orders.  Equation  (9)  is  good  only  to  first  order 
in  n ] . ^ 

Special  cases  are  evaluated  below  for  comparison  with  pre- 
viously published  results: 


n • 1 


_ - - - 

H(pi,  P2;  Pi,  P2>  ■ e 12 


Thin  is  a well  known  result, 
n - 2 


36-38 


H(pi,  P2,  P3»  P4?  Pi*  P2*  P3»  P4>  - exp  - 7 0>12  " °13  + D14 


+ I>23  - D24  + D34)  + 2C  +2C  ' 

13  24_ 


This  is  a generalization  (double  arguments)  of  the  result  of 
Ref.  39. 

Although  the  above  formulation  for  the  generalized  spherical 
wave  mutual  coherence  function  is  strictly  valid  only  when  p is 
normally  distributed,  this  is  not  an  overly  restrictive  require- 
ment. From  physical  reasoning  and  theoretical  considerations, 

it  can  be  shown  that  both  the  log  amplitude  and  the  phase  are 

29,40 

normally  distributed  to  first  order  in  nj.  ’ This  has  been 

41-49 

confirmed  by  experiment.  Also,  there  is  experimental 

evidence  to  suggest  that  x is  also  normally  distributed  or 

nearly  so  for  situations  of  practical  importance  under  con- 

29  40-47 

ditions  of  multiple  scattering.  ’ Since  the  larger  scale 

sizes  are  more  important  with  respect  to  phase  effects,  it  seems 
reasonable  that  multiple  scattering  should  be  less  important 


69 


and  that  the  first  order  theory  for  phase  effects  should  be 

valid  even  when  the  log-amplitude  fluctuations  are  highly 

saturated.  Measurements  ’ of  angle  of  arrival  fluctuations 

tend  to  confirm  this  conclusion  and  it  is  generally  relieved 

29,40 

that  <j>  i'  indeed  normally  distributed  with  <<J>>  = 0. 

Consequently,  (7)  should  be  valid  in  general  for  appropriate 

29,40 

ranges  of  log-amplitude  variance  and  valid  without 

restriction  in  special  cases  such  as  a large  Fresnel  number 

50 

transmitter  where  the  phase  term  is  dominant. 
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